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We reexamine the dimerization, the charge and the spin gaps of a half-filled Peierls- 
Hubbard chain by means of the incremental expansion technique. Our numerical findings 
are in significant quantitative conflict with recently obtained results by M. Sugiura and 
Y. Suzumura [J. Phys. Soc. Jpn. 71 (2002) 697] based on a bosonization and a renor- 
malization group method, especially with respect to the charge gap. Their approach 
seems to be valid only in the weakly correlated case. 
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1. Introduction 

The interplay of electron-electron (el-el) and electron-lattice (el-ph) interactions in the context 
of various phase transitions (instabilities) is one of the most challenging many-body problems in 
solid state physics. In particular, in the context of the longstanding debate how to describe the- 
oretically conducting polymers like trans polyacetylene and its remarkable solitonic properties, 1,2 ) 
the Peierls-Hubbard model (PHM) at half filling, considered as a minimal realistic model, became 

very popular. 2 ~ 14 ) Within this frame in the adiabatic approximation for the lattice the behavior 
of the dimerization amplitude as a function of the two remaining model parameters, the onsite 
repulsion U and the strength of the el-ph interaction is one of the central problems. To the best of 
our knowledge it was first addressed in 1981 by P. Horsch 3 ) who predicted on qualitative grounds a 
maximal dimerization at some intermediate U. This was confirmed later on and made more clear by 
approximate solutions obtained by Baeriswyl and Maki 5 ) and Kuprievich 11 - 1 within the Gutzwiller 
approximation and the geminal techniques, respectively. According to them the maximal dimer- 
ization occurs for weak el-ph interaction at U ~ 4to, where to denotes the nearest neighbor transfer 
integral. Unfortunately, the region of interest is already at the border line of validity of those 
approaches. 

The general nonmonotonous behavior of the dimerization amplitude (at weak el-ph coupling) 
as a function of U was confirmed by Hayden and Soos 6 ) and Waas et al?> who extrapolate results of 
exact diagonalization studies (EDS) obtained for finite even-membered rings with various boundary 
conditions. In the strongly correlated limit U S> to the PHM can be reduced to the spin-Peierls 
model (SPM) for which a large body of literature exists (see for instance Ref. 15 ) and references 
therein). The general case of arbitrary U and intermediate as well as strong el-ph interactions has 
been covered by the present authors within the solitonic 12 ) as well as the incremental expansion 
method (IEM) 14 ) both based on self-consistent EDS. 

In the context of the polyacetylene problem the charge gap has been addressed among others 

by Rosciszewski and Oles 10 - 1 in the framework of the extended PHM. There is now a general 
consensus that the Peierls mechanism, i.e. the off-diagonal el-ph interaction, enhances the charge 

gap compared with that of the lD-Hubbard-Mott insulator. 7 ~ 9 ) The latter problem has been solved 

exactly by Lieb and Wu. 16 ) 
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Recently, the weak and intermediate region U < 4io has been readdressed by Sugiura and 
Suzumura 17 ) who employ the bosonization technique and the renormalization group method. Beside 
the dimerization they consider also the U dependence of the charge and the spin gaps. It is the aim 
of the present paper to show, that these calculations are in some conflict with results obtained for the 
original PHM applying the incremental technique supplemented by density matrix renormalization 

group (DMRG) calculations. 13 ) 

2. Models and Methods 

Let us consider even-membered chains with periodic boundary conditions modeled by the 
Peierls-Hubbard Hamiltonian supplemented with a static harmonic lattice term 

H = -t it i +1 [cl jS Ci +ljS + h.c] + -Un^ s ni- S 

i,s 

i 

where the bond alternation Vi has been introduced 

Vi = u i+ i-Ui, (2) 

here Ui = ±uq is the displacement of the ith lattice site in the dimerized state under consideration 
and K is the spring constant. We linearize the bond-length dependent transfer integral i^j+i 

t iji+1 = t - ryvi. (3) 

The strength of the off-diagonal el-ph interaction can be expressed by the dimensionless coupling 
constant g 



g = j/VKt < 1. (4) 
It is convenient to measure the dimerization amplitude vq = 2uq also in dimensionless units 11 ) 

d = vq^kJTq < 1. (5) 

Values of g ~ 0.4 to 0.5 are typical for conducting polymers. For the convenience of the reader we 
note that our quantities g and d are simply related to the notations A = Xssh introduced by Su, 
Schrieffer and Heeger 1 ) as well as used by Sugiura and Suzumura 17 ) 

A = 4g 2 /7r, t d = dgt . (6) 

Thus their two values of A =0.25, 0.42 correspond to our g=0. 443113 and 0.574341 (used in the 
calculation shown in Figs. 1-4), respectively. 

Let us start with the calculation of the dimerization amplitude d. Naturally in the present 
homogeneous dimerized case it can be obtained from the direct minimization of the total energy 
per site with respect to the single variable d at fixed U and g, in other words there is no need 
to solve the general self-consistency equations presented in Ref. 12 ) According to the incremental 
expansion technique the ground state energy per site in the n th order reads 14 ) 

e n (U,g,d) = 1 [E 2n+2>2n+2 (U,g,d) - E 2n>2n (U, g,d)] , (7) 

where E 2riy2n (U, g, d) is the total energy of a chain segment (finite open chain) with 2n electrons and 
2n sites. The dimerization amplitude in the infinite chain limit (N — > oo, where N = 2n + 2 is the 
maximal number of interacting sites in the n th incremental expansion order) has been extrapolated 
using the empirical law (see also 12 ) ) 

d(U,g,N) = (kofrg) + (^2 ai (U,g)/N^j exp(-£([/, g)N). (8) 
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Notice that £ is close to the coherence length (half-width of a soliton as obtained for instance by 
the solitonic method 12 ) explicitly). 

As the next step the gaps have been calculated for the extrapolated dimerization d^. In the 
n order of the incremental expansion the charge gap (for definition see Ref. 16 )) reads 

^£(n) = ^+ _ ^ _ (n + 1) [E2n+3,2n+2 ~ £?2n+2,2n+2] ~ U [^2n+l,2n — ^2n,2n] 
— (n + 1) [-E2n+2,2n+2 ~ ^2n+l,2n+2] ~ Tl [E2 n ,2n ~ ^2n-l,2n] 
= (n + 1) [E 2n +3,2n+2 — 2-E , 2n+2,2n+2 + ^2n+l,2n+2] 

— n [E2n+l,2n ~ 2£'2n,2n + ^2n-l,2n] > (9) 

where / u ± denotes the chemical potential for adding an electron or a hole, respectively. 

The spin gap A£ a can be calculated using an analogous expressions as for fi + (i.e. the first line 
of Eq. (9)) 

A£(«) = (n + 1) (-E' 2rt+ 2 l 2n+2 ~~ -E , fn+2,2n+2) ~~ n {^2n,2n ~ ^2n,2n) ) (10) 

where Et(s) ^ s the total energy of the lowest state in the = 1(5^ = 0) subspace. These states 
correspond usually to the triplet and singlet states, respectively. For a more detailed derivation of 
Eqs. (9,10) and technical details of the incremental technique the reader is referred to Refs. 14 ' 18 ) 

Notice that in the expression for the charge gap (Eq. (9)) the corresponding multiplet notation 
(doublet or singlet) has been omitted for the sake of shortness. Let's note also that in Ref. 17 ) the 
charge and spin gaps have been denoted by A p and A a , respectively. 

The Peierls contribution to the total charge gap as a function of U has been determined using 

relent ion. 

5£jr\U,g) = A£p\u,g)-A£<?>(U,0), (11) 

where A£ p n \ll, 0) denotes the charge gap of the undimerized chain calculated in the n th incremental 
order. The quantity 5£ p n \u,g) has then extrapolated by a polynomial fit with the leading term 
oc iV~ 2 to the infinite chain limit 8£ p (U,g). Finally, 5£ p {U,g) has been added to the charge gap 
of the equidistant undimerized chain (sometimes called also Mott-Hubbard gap) given by the well- 
known Lieb-Wu solution 16 ) 

oo 

A£ p (U, 0) = U - 4 - 8 J^(-l) n + (nU/2) 2 ) - nU/2\ . (12) 

n=l 

This way we circumvent some difficulties with the size dependence of the total charge gap which 
stem mainly from the slow size dependence of the Mott-Hubbard gap at small U. 19 ^ 

3. Results and discussion 

Let us start with the discussion of our numerical results obtained up to 6 th order in the incre- 
mental expansion using the method of exact diagonalization for the Peierls-Hubbard Hamiltonian 
of the chain segments. All results have been extrapolated as mentioned in the previous section. 
Dimerization amplitudes t& are shown in Fig. 1. From the comparison with Fig. 4 of Ref. 17 ) one 
realizes reasonable agreement in the weak coupling region below 1.5 U/to but a nearly 100 % larger 
dimerization amplitude at U/to=4 in our calculation. This is mainly caused by the shift of the 
maximum position U max [g] to higher U- values and to a more pronounced asymmetric shape of the 

whole curve (see Fig. 1, right panel). Notice that in Fig. 4 of Ref. 17 ) there is essentially no shift of 
the maximum dimerization towards smaller [/-values with increasing el-ph interaction A as it should 
be. In the weak coupling limit it has to tend to U max (0) w 3to whereas U max tends to zero for 
9 9c ~ 0.69. This is in contrast with the statement of U max =2.2 ±0.1 for 0.2< A < 0.5 made in 
Ref. 17 ) Furthermore, the difference 5d = d(U max ) — d(0) decreases with increasing g contrarily with 
the results shown in Fig. 4 of Ref. 17 ) where it increases by a factor of two. The same problem is also 
evident from the inspection of Fig. 3 of Ref. 17 ) For A approaching the critical value A c «0. 606189 
mentioned above, the td(U = 0) curve should intersect all t^-curves at finite U, because for all 
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cl-ph constants A > A c the maximal dimerization is already reached in the uncorrelated limit U =0. 
(Strictly speaking, their figure Fig. 3 shows only A <0.6, i.e. slightly below A c . Anyhow, no visible 
change in the slope of the id(A)-curves with finite U approaching the closely lying critical point A c 
can be detected.) 

The charge and the spin gaps are shown in Figs. 2,3. One readily realizes (i) that 5£ p (U,g) is 
always positive, i.e. the Peierls mechanism enhances the total charge gap and (ii) it follows closely 
the spin gap (see Fig. 3). Roughly speaking, there are three regions: the weakly coupled one 
U/to <1.5, the intermediately coupled and the strongly coupled ones. In the first one correlation 
effects can be neglected and the charge gap being almost coinciding with the spin gap is dominated 
by the el-ph interaction. In the second intermediate region the Peierls mechanism yields a non- 
negligible contribution to the total charge gap. This is of interest for conducting polymers, where 
U ~ U max ~ 3io- Here this Peierls contribution amounts about 1/3 of the total charge gap, 
whereas in the strongly coupled spin-Peierls region it becomes increasingly small reflecting this 
way the asymptotical decoupling of spin and charge degrees of freedom. Quite interestingly, the 
dimerization, the Peierls-contribution to the charge gap as well as the spin gap exhibit a very 
similar dependence on U. From this observation one might suggest that all those three quantities 
are mainly described by the spin degrees of freedom (see also Ref. 12 ) where an effective spin-Peierls 
model has been proposed to describe the dimerization d(U,g) even analytically). 

From our calculations it becomes clear that there is no visible hump in the charge gap curve 
near the region where the dimerization amplitude td(d) and the spin gap reach their maximum. 

More strikingly, in Ref. 17 ** the calculated charge gap at U = 4io is even smaller than the Lieb-Wu 
predictions (see Eq. (12)). This would imply a negative impact of the Peierls mechanism in sharp 
contrast with our findings (see Fig. 3). 

Finally, below we report also on some incremental expansion results where due to the ap- 
plication of the DMRG-method to open chain segments rather high orders could be considered. 
Furthermore, charge gaps obtained from total energy differences of long open chains (denoted as 
standard open chain results) will be provided, where the well-known definition of the charge gap 16 ) 

A£ p (U, g, N) = E(+e, U, g) - 2E(0, U, g) + E(-e, U, g, N) (13) 

has been used. Here E(0,U, g, N) and E(±e,U, g, N) denotes the total energy of the neutral and 
charged open chain, respectively, at fixed d. 

For instance, for g = 0.4 and U/to = 4 in the 31 st incremental order one arrives at d =0.10591 
(0.10590). With that dimerization we obtained in the 63 rd order the charge gap A£ p /to=lA992 
(1.5195) whereas the standard open chain result for N=128 sites is A£ p /to= 1.5077. Since the charge 
gap of the undimerized chain calculated from Eq. (12) equals to 1.2867 one realizes a significant 
enhancement of the charge gap due to the finite dimerization. Here and below, the numbers in 
brackets denote the results obtained by the extrapolation of the incremental data calculated up 
to the 6 th order. At the same incremental orders and open chain lengths as mentioned above, for 
g = 0.45 (which is very close to A = 0.25 considered in Ref. 17 )), and for U/to=2 we found d = 0.13575 
(0.13650) and A£ p /£ =0.5208 (0.5127) whereas the standard open chain result is A£ p /t =0.5374 
(compare with 0.1728 from Eq. (12)). For g = 0.5 and the same U/to=2 as above we have now 
d = 0.18578 (0.18580) and A£ p /to=0.6696 (0.6764) whereas the standard open chain result is 
A£" p /t =0.6828. For the same g = 0.5 but U/t =7, one arrives at d=0. 11977 and A£ p /t =3.9734 
(3.9925). The standard open chain result is A£ p /to=3.9777 whereas Eq. (12) yields 3.7720, only. 

Thus the open chain results (Eq. (13)) provide an upper bound whereas the incremental tech- 
nique (Eq. (9)) results in a lower bound for the charge gap. The incremental technique yields 
excellent results for the gaps in the intermediate and strongly correlated regime whereas in the 
weakly correlated regime slight deviations occur, if one restricts oneself to sixth order dictated by 
the computational facilities available for exact diagonalization at present. Due to an exponential 
size dependence in the dimerization amplitude (see Eq. (8)) it gives excellent results even here. 
Anyhow, if the incremental technique is combined with the DMRG technique, the higher orders 
necessary for the weakly correlated cases can be calculated with high enough accuracy. In other 
words, the incremental expansion technique combined with DMRG method provides a very pow- 
erful tool to evaluate with high accuracy the dimerization amplitude, the charge gap as well as the 
spin gap. 
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To conclude, from the critical comparison of our discrete calculations with the bosonization 
(continuum) approach proposed by Sugiura and Suzumura 17 ) one has the impression that the 
proposed formalism is useful only in the weak correlation region < U /to <l-5. It might be inter- 
esting to elucidate to what extent such an approach might be improved for instance in the spirit 
of the finite-band approach proposed by Gammel 20 ) to improve the shortcomings of the continuum 
models 21 ' 22 -* compared with the noninteracting discrete SSH-model. 1 ) Alternatively, an improved 

and extended effective spin-Peierls model based on the description proposed in Ref. 12 ) would be 
noteworthy, too, at least in the intermediately strong correlated case where the maximum of the 
dimerization occurs for weak el-ph interaction. 
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FIGURE CAPTIONS 



Fig.l 

Alternating part = gdto (dimerization) of the transfer integral t vs. on-site (Hubbard) repulsion 
U for the two el-ph coupling strengths considered in Ref. 17 -* (left panel). The same as in the left 
panel for a larger scale of U and slightly changed el-ph coupling constants g (right panel). 

Fig. 2 

The Peierls contribution to the total charge gap (see Eqs. (8,10) and the spin gap (Eq. (9) vs. on-site 
repulsion U. For comparison the charge gap of the equidistant infinite chain, i.e. the Mott-Hubbard 
gap (full line, no symbols) is shown too. 

Fig. 3 

The Peierls contribution to the total charge gap and the spin gap vs. on-site repulsion U. 
Fig. 4 

Comparison of extrapolated and 6 th -ovder incremental calculations. 
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